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ABSTRACT 



Wc measure the apparent shapes for a sample of 62 blue compact dwarf 
galaxies (BCDs), and compare them with the apparent shapes for a sample 
of 80 dwarf elliptical galaxies (dEs). The BCDs are flatter, on average, than 
the dEs, but the difference is only marginally significant. We then use both 
non-parametric and parametric techniques to determine possible distributions of 
intrinsic shapes for the BCDs. The hypothesis that BCDs are oblate spheroids 
can be ruled out with a high confidence level (> 99%), but the hypothesis 
that they are prolate spheroids cannot be excluded. The apparent shapes of 
BCDs are totally consistent with the hypothesis that they are triaxial ellipsoids. 
If the intrinsic axis ratios, f3 and 7, are distributed according to a Gaussian 
with means and 70 and standard deviation a, we find the best-fitting 
distribution for BCDs has (/5o, 70, cr) = (0.66, 0.55, 0.16), while that for dEs has 
(/3o,7o,cr) = (0.85,0.64,0.24). Our results are consistent with the hypothesis 
that BCDs have a close evolutionary relation with dEs. 



Subject headings: galaxies: blue compact dwarf - galaxies: dwarf elliptical - 
galaxies: structure 
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1. Introduction 

Understanding the formation and evolution of stellar systems like galaxies begins from 
the comparison of various observables. Dwarf galaxies are the most numerous type of galaxy 
in the universe. Nearly 80% of the known Local Group galaxies are dwarfs, and the number 
density of dwarf galaxies might be several dozen times that of giant spirals and ellipticals 
(Staveley-Smith, Davies, & Kinman 1992). The most common type of dwarf galaxy is the 
dwarf elliptical (hereafter dE). They have very smooth surface brightness profiles, which 
are typically nearly exponential (Sandage & Binggeh 1984; Caldwell & Bothun 1987). The 
dEs, like giant ellipticals, are found predominantly in groups and clusters of galaxies. Most 
dEs are symmetric and apparently relaxed with very uniform colors indicating that they are 
dominated by old stars (James 1994). The best-studied dwarf ellipticals are those in the 
Virgo cluster, cataloged by Sandage & Binggeli (1984) and Binggeli, Sandage, & Tammann 
(1985). 

A smaller group of dwarf galaxies are classified as "blue compact dwarfs" (hereafter 
BCDs). In contrast to dEs, BCDs are undergoing intense bursts of star formation which 
give birth to a large number of massive stars, as evidenced by their blue UBV colors. These 
very young hot stars ionize the interstellar medium, producing high-excitation super-giant 
HII regions (Thuan 1991; Thuan, Izotov, & Lipovetsky 1995). As a result, the optical 
spectra of BCDs show strong narrow emission lines. Their ultraviolet spectra show a steeply 
rising continuum toward the blue, which is characteristic of OB stars (du Puy 1970; Searle 
& Sargent 1972). 

In addition to these spectroscopic and spcctrophotometric observations, another 
important tool for the comparison of different populations of dwarf galaxies is provided 
by morphological studies - for instance, the study of intrinsic shapes based on measured 
ellipticity distributions. There have been several trials which have explored this point: 
Caldwell (1983), Ichikawa, Wakamatsu, & Okamura (1986), Ichikawa (1989), Ferguson & 
Sandage (1989), and Ryden & Tcrndrup (1994) for dEs and Staveley-Smith et al. (1992) for 
BCDs. However, due to different observational techniques, different methods of measuring 
ellipticities, and different algorithms for deriving intrinsic shapes, direct comparison of the 
results might cause systematic errors. In addition, some previous studies suffered from 
additional systematic errors in eye-estimated flattening data (Fasano & Vio 1991; Ferguson 
& Binggeli 1994). Among the above authors, only the analysis of Ryden & Terndrup (1994) 
was based on CCD observations of a large data set. The analysis of Staveley-Smith et al. 
(1992) was also based on CCD observations; however, they determined the intrinsic shapes 
of BCDs under the (probably erroneous) assumption that all BCDs are oblate spheroids. 
In addition, the total size of their sample is too small to draw any strong conclusion. For 
studies of intrinsic shapes, a large homogeneous data set is essential. 

In this paper, we measure the ellipticities of 62 BCDs, and compare the ellipticity 
distribution with that of a sample of 80 dEs (Ryden & Terndrup 1994; Ryden et al. 1997). 
The ellipticity distribution of BCDs proves to be marginally different from that of dEs. We 
then determine the distribution of intrinsic shapes for the BCDs, under various assumptions. 
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First, we assume that BCDs are either all oblate or all prolate, and use non-paramctric 
analysis to find the best-fitting distribution of intrinsic shapes. Next, we assume that BCDs 
are triaxial, and use parametric analysis to find permissible distributions of intrinsic shapes. 
Prom these analyses, we find that the shapes of BCDs are well described either by triaxial 
ellipsoids or prolate spheroids; the hypothesis that BCDs are all oblate is ruled out. Similar 
results are obtained for dEs. If the intrinsic axis ratios, (3 and 7, are distributed according 
to a Gaussian with means /3o and 70 and a common standard deviation of a, we find the 
best-fitting intrinsic axis ratio distribution of BCDs has (/5o, 70, cr) = (0.66,0.55,0.16) with 
the KS probability of Pks = 0.99, while that of dEs has (/3o, 7o, cr) = (0.85, 0.64, 0.24) with 
Pks = 0.97. Given the similarity of shape distributions, our results are consistent with the 
hypothesis that BCDs have a close evolutionary relation with dEs. 



2. Observations 

To measure the apparent axis ratios, we obtained images of 62 BCDs with the 1 m 
telescope at Mount Stromlo & Siding Spring Observatory. Eight observing sessions were 
carried out from 1993 April to 1994 March, employing a 1024 x 1024 Tektronics chip, 
which provides a field of view of 10' x 10' with a pixel size of 0".57pixeP^. The exposure 
time varied from 300 s to 600 s depending on the brightnesses of individual galaxies. We 
attempted to obtain the 4 (Bessell 1990) BVRI band images of individual BCDs. However, 
due to limited observing time, we obtained full 4-band images for only 35 out of 62 galaxies. 
For the other galaxies we obtained partial band images: 3 bands for 7 galaxies, 2 bands for 
another 7 galaxies, and only a single band image for the remaining 13 galaxies. Therefore, 
the total number of images is 188. The observed bands of each image are listed in Table 1. 

We selected our sample of BCDs from various sources. We choose nearly half of our 
sample galaxies from the list of Thuan & Martin (1981); and the rest are from Maza et al. 
(1991), Bergvall & Olofsson (1986), Kunth & Sargent (1986), Macalpine & Williams (1981), 
Fairall (1977), Terlevich et al. (1991), Gondhaeker et al. (1984), Wamsteker et al. (1985), 
and Acker, Stenholm, &: Verson (1991). The galaxies in these sources are classified as BCDs 
according to 3 major criteria: 

1. low luminosity, with Mb j:^ — 18, 

2. strong, sharp, narrow emission lines superimposed on a blue continuum, 

3. very compact, with diameters of a few kpc. 

To meet the necessities of observation and data analysis, we impose several additional 
criteria: 

4. relatively small recession velocities [v ^ 4, 000 km s~^ ), 

5. relatively bright apparent magnitudes {B < 16), 

6. relatively large angular sizes 20"), 

7. unperturbed outer isophotes. 

We include conditions 4 and 5 to select galaxies bright enough to obtain high S/N images 
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with the available telescope. Condition 6 is included because it is difficult to measure 
ellipticities precisely for galaxies of small angular size. According to the classification scheme 
of Loose & Thuan (1986), BDCs are classified into four types based on the regularity of 
the isophotal shapes of both the high surface brightness inner region, where star formation 
has recently occurred, and the low surface brightness outer region, which contains only an 
older stellar population. The first and most common type, comprising ~ 70% of BCDs, 
is 'iE', which combines irregular inner isophotes with accurately elliptical outer isophotes. 
The 'nE' type has elliptical isophotes in both the inner and outer regions. The 'il' type has 
irregular inner isophotes and outer isophotes that deviate measureably from perfect ellipses. 
Finally, 'iO' galaxies do not show faint outer structures at all; they appear to contain only 
starburst regions, and lack an underlying old stellar population. In our sample, we include 
all BCDs as long as their outer isophotes can be approximated as ellipses; our sample thus 
contains only iE, nE, and il galaxies. However, we note that our sample represents the 
whole BCD population well because BCDs of type iO are very rare. 

All steps of the data reduction and analysis were performed using a standard CCD 
reduction process with the Image Reduction and Analysis Facilities (IRAF). Raw images 
were first corrected for zero level bias and overscan pattern. Each image was then divided 
by a fiat field obtained from the median of 5 - 10 twilight and dome flat images. During 
this process, pixels in the fiat field image which are more than 1.5 times the median 
or less than 0.5 times the median are replaced with the median value. In addition, we 
combined images of individual galaxies if multiple images taken at the same bands are 
available. By combining images, we could not only correct for cosmic rays, but also increase 
signal-to-noise ratios. In the next stage, the sky brightness was subtracted from individual 
images. For sky subtraction of most BCD images, we used separate sky images taken just 
after the observation of each galaxy. However, if no sky image is available, we used the sky 
value in the vicinity of a galaxy within the same image. Finally, the processed images are 
trimmed centered at the galaxies. Since most BCDs in our sample have angular sizes less 
than 2', the typical trimmed images have sizes of ^ 500 X 500 pixeP. If stars are located on 
or near the galaxies, we remove them by replacing them with the median value of the area 
near the stars. 



3. Determination of Apparent Axis Ratios 

To determine intrinsic axis ratio distribution, one should start by measuring the 
apparent axis ratio, g = 1 — e = 6/a, of galaxies, as projected onto the sky. Here a and h are 
the semimajor and semiminor axes of a projected galaxy, and e is the ellipticity. From the 
obtained image of each BCD, we measure the ellipticities as a function of semimajor axis 
from the center, using the Space Telescope Science Data Analysis System (STSDAS) routine 
'isophotc'. In the routine, each galaxy was modeled as a nested series of ellipses, yielding 
the semimajor axis, a, position angle, 6'pos, surface brightness, S(a), and the positions of 
each eUipse's center. 
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The cllipticitics of the isophotes of many stellar systems vary as a function of semimajor 
axis. For BCDs wc have the same problem of varying ellipticities, especially in the inner 
regions of BCDs. The main reason that makes the isophotes of a BCD deviate from 
well-aligned concentric ellipses is the existence of inner substructures, mainly bright HII 
regions (Thuan & Martin 1981; Gordon & Gottesman 1981). In Figure 1 we present, as 
an example, the B-band isophotes of ESO 495-G21, upon which are superimposed the 
best-fitting ellipses. To better show the poor ellipse fitting in the inner region of the galaxy, 
we expand the part enclosed by a box and present it in the upper right panel, in which the 
isophotes and corresponding best-fitting ellipses at a = 14.5 and 16.5 are shown. We also 
present the ellipticity as a function of the semimajor axis in the lower panel. (Note that 
all length scales are in units of the pixel width.) The isophotes are distorted by a bright 
HII region located to the right of the galaxy's center. As a result, the obtained best-fitting 
ellipses are more elongated than those which would be obtained without the HII region, and 
the position angle is tilted toward the direction of the HII region. Beyond a certain point 
(a ~ 30 pixel for ESO 495-G21) there no longer exists any substructure, and the values of 
the ellipticity become stable as shown in the figure. Note that the concentric circles centered 
at {x, y) ~ (265, 220) and (x, y) ~ (262, 189) are foreground stars, before they are removed 
for ellipse fitting. The ellipse-fitting process fails as the surface brightness of a galaxy 
approaches background value; for the case of ESO 495-G21 this point corresponds to a ~ 95. 
Therefore, as a representative axis ratio of each BCD we determine the intensity-weighted 
axis ratio averaged over the intermediate region where stable ellipticities can be obtained 
with successful ellipse-fitting processes. The intensity-weighted mean axis ratio is computed 

by 



- !q{a)dL 

Q = ; dL = 2nqa 



^ ^ 1 d\nq{a) 



2 dlna 



J:{a)da, (1) 



following the method by which Rydcn & Terndrup (1994) computed the axis ratios of dEs. 
If more than a single band images are available, we take the mean value of the ellipticities 
determined in individual bands as a representative ellipticity. The finally determined 
ellipticities of our sample galaxies are listed in Table 1. We estimate the errors of e = 1 — ^ 
by computing the variance of e(a) within the range of semimajor axis where ellipticities are 
measured. 

A BCD with lots of star formation might have a B image different from its / image, 
for instance; the B image emphasizes star formation regions, and hence is much clumpier. 
Thus, the measured eUipticities at B might be systematically different from the eUipticities 
in other bands. To test this possibihty, in Figure 2 we compare the eUipticities measured 

in B band, e^, to those measured in I band, ej. Ellipticities do not vary significantly from 
one band to another. The values of ellipticities are consistent regardless of the observed 
bands because the determination of e depends strongly on the ellipticities beyond the inner 
clumpy regions. 

We also check the consistency of our measurements by measuring ellipticities in a 
different way. In the second method, we measure the ellipticity 624 at the point where 
the surface brightness is Eg = 24 mag arcsec"^. For these determinations, we carried 
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out S-band surface photometry of arbitrarily selected 16 sample BCDs by using elliptical 
isophote fitting method with the STSDAS routine 'isophote'. Photometric calibrations 
were carried out by using ~ 50 standard stars in the list of Graham (1982) and Landolt 
(1992). In Figure 3, we plot the correlation between e and €24. The ellipticities determined 
in both ways are consistent each other. In Figure 3, the error estimates of 624 are derived 
by computing the standard deviation of isophotes from their corresponding best-fitting 

ellipses, i.e., a^.,^ = [E(xisophotc - Xcmpsc)^/npoint]^^^ Here Xisophotc and XdUpsc represent the 
positions of isophotes and ellipses on CCD pixel coordinates, and npoint is the number of 
pixels on which individual curves (both isophotes and ellipses) are drawn. 

Since we use the same method for determining e as Ryden & Terndrup (1994), we 
can directly compare our results for the apparent shapes of BCDs to their results for the 
apparent shapes of dEs in the Virgo Cluster. To the data obtained for 70 dEs by Ryden 
& Terndrup (1994), we add data for an additional 10 dEs observed in Virgo by Ryden et 
al. (1997). The cumulative distribution function of g = 1 — e for the 62 BCDs is shown as 
the solid step function in Figure 4. The cumulative distribution function of q for the 80 
Virgo dEs is shown as the dotted step function in the same Figure. A Kolmogorov-Smirnov 
test comparing the two distributions shows a marginal difference, with Pks = 0.057. The 
main difference between the apparent axis ratios of the two distributions is that BCDs are 
slightly flatter, on average. For the sample of BCDs, the mean and standard deviation of q 
are 0.671 ± 0.150; for the sample of dEs, they are 0.719 ± 0.156. A Mann- Whitney U-test 
comparing the medians of the two data sets yields a probability Pmw = 0.068. Thus, we can 
state from the Mann- Whitney test that the two distributions differ at the 93.2% confidence 
level. 



4. Intrinsic Shape Determination 

The apparent axis ratio distributions of stellar systems provide important constraints 
about the 3-dimcnsional shapes of the systems (Fasano & Vio 1991; Lambas, Maddox, 
& Loveday 1992; Tremblay & Merritt 1995; Ryden 1992, 1996). For the determination 
of intrinsic shapes, two methods are often applied. Both methods assume that the 
galaxies which we observe are randomly oriented with respect to us. The first method for 
determining intrinsic shapes assumes that the galaxies are intrinsically axisymmetric, and 
that they are either all oblate or all prolate. With these assumptions, one can perform a 
unique mathematical inversion to go from the distribution f(q) of apparent shapes to the 
distribution 7(7) of intrinsic shapes. The second method assumes that galaxies are triaxial; 
in this case one no longer has a unique solution to the problem. However, if one assumes 
a functional form (e.g., Gaussian) for the distribution of intrinsic axis ratios, one can vary 
the parameters of the function until one finds the model which best fits the observed 
distribution of apparent axis ratios. Parametric fits are useful because they enable us to 
see how statistics such as the and KS scores vary as the parameters are changed. We 
estimate the intrinsic shapes of BCDs by using both the non-parametric and the parametric 
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methods. 



4.1. Non-pctrametric Method 

Using the non-parametric method, we reject or accept, at a known confidence level, the 
null hypothesis that the intrinsic shapes of the galaxies in a sample are randomly oriented 
oblate spheroids. Moreover, if the null hypothesis is not rejected, the non-parametric 
method gives an estimate for the distribution function fo of the galaxies' intrinsic axis ratios, 
7 < 1. Similarly, we can reject or accept the hypothesis that they are randomly oriented 
prolate spheroids. To accomplish this, we make a non-parametric kernel estimate f{q) of 
the distribution of the apparent axis ratio q. We then mathematically invert f{q) to find 
/o(7) and /p(7), the estimated distributions of the intrinsic axis ratio 7 for a population of 
oblate and prolate spheroids, respectively. Confidence intervals are placed on the estimators 

fo, and /p by performing repeated bootstrap resampling of the observed data and 
creating new estimates from each sample. The spread in the bootstrap estimates of / at a 
given value of q and of fo and /p at a given value of 7 provides confidence intervals for the 
non-parametric estimates of these functions. Once the estimators fo and /p are determined, 
one can exclude the oblate or prolate hypothesis at a given confidence level if the upper 
confidence limit drops below zero for any value of the intrinsic axis ratio 7. For the details 
of the non-parametric kernel estimators, see Ryden (1996) and Tremblay & Merritt (1995). 

In the upper panel of Figure 5a, we present the non-parametric kernel estimate of 
the distribution of the apparent axis ratios for our sample of 62 BCDs. In the following 
panels, we also present the distributions of intrinsic axis ratios assuming that the BCDs 
are oblate (middle) and prolate (bottom) . The solid line in each panel is the best estimate, 
the dashed lines are the 80% confidence band, and the dotted lines are the 98% confidence 
band. (That is, 10% of the bootstrap resamplings lie above the 80% confidence band, and 
10% lie below it. Similarly, 1% of the resamplings lie above the 98% confidence band, and 
1% he below it.) The hypothesis that BCDs are randomly oriented oblate spheroids is ruled 
out at the 99% (one-sided) confidence level. Due to the lack of nearly circular BCDs (those 
with q ~ 1), the 98% confidence band for fo dips far below zero at 7 > 0.9. On the other 
hand, the prolate hypothesis cannot be ruled out at the 99% (one-sided) confidence level. 
The best-fitting prolate distribution (the solid line in the bottom panel of Figure 5a) gives 
a mean and standard deviation for the intrinsic axis ratio 7 of 0.58 ± 0.13. 

For the comparison of intrinsic shape of BCDs with that of dEs we also present /(g), 
/o, and /p for the sample of 80 Virgo dEs, in Figure 5b. We find that the oblate hypothesis 
cannot be ruled out at the 99% (one-sided) confidence level for this sample, while the 
prolate hypothesis is perfectly acceptable. The best-fitting prolate distribution gives a mean 
and standard deviation for 7 of 0.63 ± 0.15. Ryden (1996), analyzing a set of 170 Virgo 
dEs observed by Binggeh & Cameron (1993), found that for this larger sample, the oblate 
hypothesis can be ruled out at the 99% confidence level. By contrast, Ichikawa (1989), 
analyzing a sample of 98 Virgo dEs, alleges that their apparent shapes are consistent with 
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their being a population of oblate spheroids. Instead of doing a non-parametric analysis, 
Icliikawa used Gaussian approximations to the correct values of /o and /p, finding that 
a Gaussian with mean and standard deviation (70, a) ~ (0.6, 0.06) is the best fit for a 
population of oblate spheroids. However, although this is the best fit, it is not a good fit. 
Although he does not provide a probability, fig 1(a) of his paper shows very poor fit; 
the best-fitting distribution overestimates the number of nearly circular dEs and severely 
underestimates the number of apparently highly flattened dEs. Thus, Ichikawa's results, 
like ours, indicate that the oblate hypothesis yields a poor fit to the observed distribution 
of apparent axis ratios. Ichikawa also presented the best-fitting distribution under prolate 
hypothesis with (70,0") ~ (0.64,0.12), which is very similar to our determination of the 
mean and standard deviation of 7 for prolate dEs. 



4.2. Parametric Method 



For the parametric determination of the intrinsic axis ratio distribution of BCDs, we 
follow the assumptions of Rydcn & Tcrndrup (1994); that galaxies arc triaxial ellipsoids 
with axis ratios 1 > /5 > 7, and that the distribution of intrinsic axis ratios follows a 
Gaussian distribution with means /3o and 70 and a common width a; i.e.. 



fiP, 7) exp 



(/3 - Pof + (7 - 70)^ 

2(72 



(2) 



With these assumptions we produce a large number (10^) of test BCDs whose intrinsic axis 
ratios are distributed according to equation (2). Then we compute their projected axis 
ratios assuming that BCDs are randomly oriented with respect to us. When an triaxial 
ellipsoid is projected with the viewing angles of 9 and 0, it appears as an ellipse with an 
apparent axis ratio of 



A+C-^{A-C)^ + B^ 



A + c + ^{A-cy + 52 



1/2 



(3) 



where 



A 
B 
C 



(cos^ + /^^ sin^ 0) cos^ 9 + ^^ sin^ 9, 
cos^sin20(l 
sin^ 



(4) 

+ COS^ 

(Binney 1985). The cumulative distribution of the projected axis ratios of the test BCDs 
is then compared with the observed one by using Kolmogorov-Smirnov (KS) tests; we 
additionally apply tests to the binned distributions. 

Figure 6a and 6b show the isoprobability contours on 4 slices through the {Po,1o, o") 
parameter space, as measured by KS and tests, respectively. When measured by KS 
tests the best-fitting distribution has parameters of {Po:lo:'^) — (0.66,0.55,0.16) with the 
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KS probability of Pks = 0.99, implying that the intrinsic shape of BCDs can be well fitted 
by a population of triaxial ellipsoids. We obtain consistent results when measured by 
tests; the best-fitting distribution has parameters of (/3o, 70, cr) = (0.77,0.51,0.16) with the 
probabihty of P^2 = 0.96. 

For the 80 dE galaxies from Ryden & Terndrup (1994) and Rydcn ct al. (1997), the 
best-fitting parametric model, as measured by a KS test, had (/5o, 70, cr) = (0.85, 0.64, 0.24), 
with Pks = 0.97. In Figure 7, we show four shces through parameter space, showing the 
isoprobability contours for the parametric dE models as measured by KS tests. When the 
best-fitting distributions for dEs arc compared to that of BCDs, one finds that the dEs are 
rounder, but the difference is not greatly significant. 

In Figure 4, we present the cumulative distributions of apparent axis ratios for BCDs 
(solid step function) and dEs (dotted step function). On top of each observed distribution we 
overlay the expected distribution from the best-fitting triaxial model, as measured by a KS 
test. For BCDs, the prediction of the best-fitting model, with {Pq, 70, cr) = (0.66, 0.55, 0.16), 
is given by the solid smooth line. For dEs, the prediction of the best-fitting model, with 
(/5o, 7o, cr) = (0.85, 0.64, 0.24), is given by the dotted smooth line. 



5. Discussion 

The results of our morphological analysis based on 62 BCDs and 80 dEs are summarized 
as follows. 

1. The apparent shapes of both BCDs and dEs can be explained by their being either 
triaxial ellipsoids or prolate spheroids. 

2. However, it is unlikely that the intrinsic shapes of either BCDs and dEs are purely 
oblate spheroids. 

3. Compared to dEs, BCDs are flatter, but the difference between the apparent axis 
ratio distributions of BCDs and dEs is marginal. 

The nature of BCDs remains ambiguous and confusing (Papaderos et al. 1996), and 
thus theories to explain their nature are very controversial. There are two major competing 
hypotheses. The first one claims that BCDs are basically a different population from dEs. 
According to this scenario BCDs are truly young systems, in which the present starburst is 
the first in the galaxy's lifetime. The other model suggests that BCDs might be dEs which 
are mainly composed of old stellar populations, and that the observed spectroscopic features 
and the spectral energy distributions arc attributed to a recent burst of star formation. 
That is, BCDs are basically dEs seen during a brief episode {t ~ 10^ yr) of violent star 
formation. 

Although strong conclusions cannot be derived just from a single comparison of 
morphology, the similarity between the intrinsic shapes of BCDs and dEs may give support 
to the second hypothesis that BCDs are just a sub-population of dEs that have undergone 
a recent burst of star formation. The second hypothesis is also supported by other 
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observations. Thuan (1983) argued that the near-infrared emission in the vast majority 
of BCDs is attributable to old K and M stars, which are the major components of dEs; 
Hunter & Gallagher (1985) also agreed that young stars are not the major contributors to 
the near-infrared flux of BCDs. In addition, many BCDs, like dEs, have surface brightness 
profiles that are well fitted by exponentials. For example, Haro 2 was classified as a BCD 
by Haro (1956) and Markarian (1967), based on its spectral features. By contrast. Loose 
& Thuan (1986) claimed that Haro 2 should be classified as a dE, based on its surface 
brightness profiles. 

The question of how to explain the difference in spectra between the two types of 
dwarf galaxy still remains. Drinkwater & Hardy (1991) stated that BCDs, after the end of 
starburst activity, fade to form dEs, as the amount of gas decreases by stripping. On the 
other hand, Silk, Wyse, & Shields (1987) proposed a totally reversed evolution sequence: 
from dEs to BCDs. According to them, hot gas, left over from the galaxy formation era, 
is trapped in a galaxy group, cools and accretes onto dEs, and provokes an intense burst 
of star formation. As a result, dEs are transformed into dEs. Although our results cannot 
clearly designate the direction of evolution, they suggest some evolutionary relation between 
the two systems. 

We would like to thank staffs at the Mt. Stromlo & Siding Spring Observatory, 
especially K. C. Freeman for supporting observations. We also thank D. M. Terndrup, 
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ellipticals. E.-C. S. has been supported by Basic Research Fund of Korea Astronomy 
Observatory. B. S. R. was supported by grant NSF AST-93-577396. 



- 12 - 



REFERENCES 

Acker, A., Stenholm, B., & Verson, P. 1991, A&AS, 87, 449 
Bergvall, N., & Olofsson, K. 1986, A&AS, 64, 469 
Bessell, M. S. 1990, PASP, 102, 1181 

Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AJ, 90, 1681 

Binggeli, B., & Cameron, L. M. 1993, A&AS, 98, 297 

Binney, J. 1985, MNRAS, 212, 767 

Caldwell, N. 1983, AJ, 88, 804 

Caldwell, N., & Bothun, G. D. 1987, AJ, 94, 1126 

Drinkwater, M., & Hardy, E. 1991, AJ, 101, 94 

du Puy, D. L. 1970, AJ, 75, 1143 

Fairall, A. P. 1977, MNRAS, 180, 391 

Fasano, G., & Vio, R. 1991, MNRAS, 249, 629 

Ferguson, H. C., & Sandage, A. 1989, ApJ, 346, L53 

Ferguson, H. C., & Binggeli, B. 1994, A&A Rev., 6, 67 

Gondhaeker, P. M., Morgas, D. H., Dopita, M., & Phillip, A. P. 1984, MNRAS, 209, 59 
Gordon, D., & Gottesman, S. T. 1981, A J, 86, 161 
Graham, J. A. 1982, PASP, 94, 244 

Haro, G. 1956, Bol. Obs. Tomantzintla y Tacubaya, 14, 8 
Hunter, D. A., & Gallagher, J. S. 1985, ApJS, 58, 533 
Ichikawa, S.-I. 1989, AJ, 97, 1600 

Ichikawa, S.-L, Wakamatsu, K.-L, & Okamura, S. 1986, ApJS, 60, 475 

James, P. A. 1994, MNRAS, 269, 176 

Kunth, D., & Sargent, W. L. W. 1986, AJ, 91, 761 

Lambas, D. S., Maddox, S. J., & Loveday, J. 1992, MNRAS, 258, 404 

Landolt, A. U. 1992, AJ, 104, 340 

Loose, H. H., & Thuan, T. X. 1986, in Star- Forming Dwarf Galaxies and Related Objects, 
ed. D. Kunth, T. X. Thuan, & T. T. Van (Paris: Editions Frontieres), 73 

Macalpine, G. M., & WiUiams, G. A. 1981, ApJS, 45, 113 

Markarian, B. E. 1967, Astrofizika, 3, 55 

Maza, J., Ruiz, M. T., Pcna, M., Gonzalez, L. E.,& Wischnjewsky, M. 1991, A&AS, 89, 389 
Papaderos, P., Loose, H.-M., Fricke, K. J., & Thuan, T. X. 1996, A&A, 314, 59 



-13- 



Ryden, B. S. 1992, ApJ, 396, 445 
Ryden, B. S. 1996, ApJ, 461, 146 
Ryden, B. S., & Terndrup, D. M. 1994, ApJ, 425, 43 

Ryden, B. S., Terndrup, D. M., Pogge, R. W., & Lauer, T. R. 1997, ApJ, submitted 

Sandage, A., & Binggeli, B. 1984, AJ, 89, 919 

Searle, L., & Sargent, W. L. 1972, ApJ, 173, 25 

Silk, J., Wyse, R. F. G., & Shields, G. A. 1987, ApJ, 322, L59 

Staveley-Smith, L., Davies, R. D., & Kinman, T. D. 1992, MNRAS, 258, 334 

Terlevich, R., Melnick, J., Masegosa, J., Moles, M., & Copetti, M. V. F. 1991, A&AS, 91, 
285 

Thuan, T. X. 1983, ApJ, 268, 667 

Thuan, T. 1991, in Proceedings of the Massive Stars in Starbursts, ed. C. Leitherer, N. 
R. Walborn, T. M. Heckman, & C. A. Norman (Cambridge: Cambridge University 
Press), 183 

Thuan, T. X., Izotov, Y. I., & Lipovetsky, V. A. 1995, ApJ, 445, 108 
Thuan, T. X., & Martin, G. E. 1981, ApJ, 247, 823 
Tremblay B., & Merritt, D. 1995, AJ, 110, 1039 

Wamstckcr, W., Pricto, A., Vitores, A., Schuster, H. E., Danks, A. C, Gonzalez, R., & 
Rodriguez, G. 1985, A&AS, 62, 255 



This manuscript was prepared with the AAS IM^jX macros v4.0. 



-14- 



TABLE 1 

The Intensity Weighted Mean Ellipticities 
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NOTE. — The intensity weighted mean ellipticities of 62 BCDs. Also hsted are the 
observed bands for each image. 
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Figure 1: Upper panels: the fi-band isophotcs of ESO 495-G21, with tlic best-fitting 
ellipses superimposed. To better show the ellipse fits for the inner isophotes, the region 
enclosed by a box is expanded and shown in the upper right panel, where the best-fitting 
ellipses at a = 14.5 and 16.5 and corresponding isophotes are presented. Note the length 

scales arc given by pixel values, which correspond to 0".57/pixel~^. Lower panel: the 
determined cllipticity of ESO 495-G21 as a function of scmimajor axis. Error estimates for 
e(a) are derived by computing the standard deviation of isophotes from their best-fitting 
ellipses. 
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Figure 2: The correlation between ellipticities determined in B and / bands. The sohd 

Une across the panel represents the identity between the two measurements, and is not a fit 
to the data. The errors in each band are estimated by computing the standard deviations 
of isophotes from the best-fitting ellipses. 
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Figure 3: The correlation between the intensity- weighted eUipticity, e, and the eUipticity 
of the Eg = 24 mag isophote, €24, for a randomly selected sample of 16 BCDs. The solid 
line across the panel represents the identity between the two measurements. The errors of 
£24 are estimated by computing the standard deviations of isophotes from the best-fitting 
ellipses, while those for e are determined from the variance of e(a). 
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Figure 4: The cumulative distribution of apparent axis ratios for 62 BCDs (solid line) 
and 80 dEs (dotted line). The step function, in each case, is the data; the smooth curve 
superimposed is the predicted distribution of apparent shapes from the best-fitting triaxial 
model. 
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Figure 5a: The non-parametric kernel estimate of the distribution of a sample of 62 BCDs 
(top panel). Also shown are the distributions of intrinsic axis ratios, assuming BCDs are all 
oblate (middle panel) and all prolate (bottom panel). The solid line in each panel is the 
best estimate, the dashed lines are the 80% confidence band, and dotted lines are the 98% 
confidence band. 
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Figure 5b: The non-parametric kernel estimate of the distribution of a sample of 80 dEs 
(top panel). The distributions of intrinsic axis ratios under all oblate (middle panel) and 
all prolate (bottom panel) are also shown. The confidence bands are same as in Figure 5a. 
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Figure 6a: Isoprobability contours, as measured by KS tests, on 4 slices through (/3o,7o,cr) 
parameter space. The fits are for the 62 galaxies in the BCD sample. Contours are drawn 
at the levels Pks = 0.01, 0.1, 0.5, and 0.9 from outside. 
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Figure 6b: Isoprobability contours, as measured by tests, on 4 slices through (/3o,7o,cr) 
parameter space. The fits are for the 62 galaxies in the BCD sample. Contours are drawn 
at the levels P^2 = 0.01, 0.1, 0.5, and 0.9 from outside. 
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Figure 7: Isoprobability contours, as measured by KS tests, on 4 slices through (/3o,7o,c) 
parameter space. The fits are for the 80 galaxies in the dE sample. Contours are drawn at 
the levels Prs = 0.01, 0.1, 0.5, and 0.9. 



